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Possible universal dynamics of a many-body system far from thermal equilibrium are explored. 
A focus is set on meta-stable non-thermal states exhibiting critical properties such as self-similarity 
and independence of the details of how the respective state has been reached. It is proposed that 
universal dynamics far from equilibrium can be tuned to exhibit a dynamical phase transition where 
these critical properties change qualitatively. This is demonstrated for the case of a superfluid two- 
component Bose gas exhibiting different types of long-lived but non-thermal critical order. Scaling 
exponents controlled by the ratio of experimentally tuneable coupling parameters offer themselves 
as natural smoking guns. The results shed light on the wealth of universal phenomena expected to 
exist in the far-from-equilibrium realm. 
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The concept of universality has been extremely suc- 
cessful in classifying and characterising equilibrium states 
of matter. For example, there are different types of or- 
der in a magnetic material separated by a second-order 
phase transition at which the relevant physical proper- 
ties become independent of the microscopic details of 
the system. This constitutes universality and allows to 
characterise an extensive range of different phenomena in 
terms of just a few classes governed by the same critical 
properties. In view of the intensifying discussion on the 
dynamics of many-body systems it is a pressing ques- 
tion whether also far away from the thermal limit the 
character of dynamical evolution can become indepen- 
dent of the microscopic details. For a closed system this 
would imply that in approaching a critical configuration 
the evolution must become independent of the particu- 
lar initial state the system has started from and critical 
slowing down in the actual time evolution is observed. 
As a consequence, different types of dynamical evolution 
could be distinguished by means of universality classes. 

In this article, we demonstrate that such universal dy- 
namics far from thermal equilibrium states is indeed pos- 
sible. We show that there is a dynamical phase transi- 
tion between different types of universal dynamics. After 
having quenched a two-component ultracold Bose gas we 
follow the ensuing evolution of the closed system and 
observe transient universal behaviour. In particular, we 
identify a tunable external parameter that determines 
the type of transient non-equilibrium order. In our setup 
non-equilibrium order is constituted by the appearance 
of spatial patterns including quasi-topological defects like 
vortices and skyrmions which are created through insta- 
bilities and which decay only on very large time scales yp . 
We find the dynamics of diluting defects to constitute a 
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FIG. 1: Illustration of the dynamical evolution, with a focus 
on the possible transition between far-from-equilibrium tran- 
sient states. The set in the bottom left corner represents an 
ensemble of out-of-equilibrium initial states. Amongst these, 
a subset subjected to a dynamical instability is marked in 
white. It is chosen such that for all values of the the ratio 
a of inter- to intra-species couplings the system belongs to 
the unstable subset. Subsequently, the time evolution of the 
closed system leads to non- and quasi-topological defect for- 
mation. The type of defects created determines the scaling 
in the low-momentum regime of the particle spectrum of the 
gas. Our results reveal three different scalings corresponding 
to three types of defect configurations, and a dynamical phase 
transition between the far-from-equilibrium universal states. 



separate dynamical critical phenomenon far from ther- 
mal equilibrium and therefore clearly distinct from the 
known equilibrium critical points [2]. Large-scale cor- 
relations in the system are universal in the sense that 
they are fixed by the type of defect, but are completely 
insensitive to the specific positions and velocities of the 
defects. Distinguishable types of defects are produced 
for different values of the tuning parameter. Hence, the 
dynamics can be tuned to a transition between different 
metastable non-equilibrium ordered states, see Fig. [l] 
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FIG. 2: Three snapshots (columns) of the evolving two-component Bose gas for three different values of the coupling ratio a 
(rows). In each case we show the spatial distribution of density imbalance as well as the momentum distribution n{k) as a 
function of the momentum modulus k of particles in the gas. Tick labels for the spectra graphs are the same as in Fig. [3] Until 
ti the dynamics is characterised by the instability and isotropisation of density fluctuations. At tNXFP = 10 ti, one observes the 
development of different power laws in the momentum distributions. They remain met ast able for a long time beyond tp = 100 ti 
and reveal different non-thermal fixed points. The change of power laws at a = 1 indicates a dynamical phase transition. 



Extending upon the presented results will allow to 
learn about the universal properties of classes of models 
in the vicinity of attractors, non-thermal fixed points and 
other critical submanifolds of the greater space of all pos- 
sible non-equilibrium configurations [SHS] . Knowing such 
universal properties, predictions for the behavior of very 
different physical systems can be obtained on the basis 
of comparatively few exemplary measurements. When 
looking at fundamental science applications, universal- 
ity classes of far-from-equilibrium criticality can link the 
dynamics, which we propose to be observable in ultra- 
cold atomic or exciton-polariton gases, with phenomena 
as different as, e.g., magnon gases in solids, quark-gluon- 
plasma dynamics in heavy-ion collisions, or reheating af- 
ter early-universe inflation. Technological applications 
would not stay away. 

The two-component Bose systems we consider are de- 
scribed by the Hamiltonian density {h = 1) 




Here m is the mass of the atoms and a = gi2/g the ra- 
tio between the inter-species coupling gi2 and the intra- 
species interaction constant g. The latter as well as the 
mass are chosen to be the same for both components and 
the sum over the field index j G {1,2} is implied. The 
considered systems allow for good experimental control 
and have been studied intensely [THIQ]. In experiment, a 
can be varied by means of a Feshbach resonance. Two dif- 
ferent ground states exist, depending on the value of the 
parameter a [Tl]|T2]. In the immiscible regime, a > 1, 
the inter-species interaction energy is greater than the 



intra-species interactions. Hence, in the ground state, the 
spatial overlap of the components is minimised through 
domain formation and spatial separation of the two com- 
ponents. On the contrary, for a < 1, the two components 
become miscible and, in the ground state, uniformly dis- 
tributed over the whole volume. We make use of a to 
change the properties of the system in the yet unexplored 
region of non-equilibrium quasi-stationary states. A fo- 
cus is set on long-lived states with non-/quasi-topological 
defects including domain walls [71 |TT] , single-species vor- 
tices, and skyrmions \i3\ in the coupled system. 

Since the dynamics we are interested in exclusively af- 
fects the low- momentum, strongly populated field modes, 
we employ the so-called classical field method which 
yields, within numerical accuracy, exact results for the 
time-evolving observables. In order to discuss the con- 
tribution of the different (quasi-) topological configura- 
tions in detail, we make use of the spin representation 
of the two-species fluid. In this representation the S^- 
component of the spin vector simply refers to the density 
imbalance between the two species. This is an essential 
observable for the detection of domain walls. Further- 
more, we make use of a hydrodynamic decomposition of 
the superfluid flow which allows for the detection of vor- 
tex contributions to the spectra. For details on numeri- 
cal parameters and decomposition methods we refer the 
reader to the supplementary material [M]. 

In Fig. § we unravel the time evolution of a two- 
component Bose gas as outlined in Fig. [T] The parame- 
ters of our simulation are chosen such that the final states 
are close to the ground state of the system. To clarify 
the type of non-equilibrium order during the intermedi- 
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FIG. 3: Decomposition of correlation functions of the two-component Bose gas in the three different non-equihbrium ordered 
states parametrised by the couphng ratio ol. Shown are the z-spin component n^, the sum of x- and y-spin components nj^ 
and the incompressible component Ui at time tNTFP • For details of the decomposition procedure see [ij . 



ate stages of the evolution, we show the imbalance be- 
tween the two components as well as the momentum 
distribution n(/c) for three characteristic times. The ini- 
tial time, ti, marks the stage of an isotropic momentum 
distribution which is overpopulated within an intermedi- 
ate range of momenta, as compared to the ensuing equi- 
librium distribution. It is marked by a strong fall-off 
at large momenta. In our driving scheme this state is 
reached in the wake of an instability. For the immis- 
cible case, a > 1, we use a modulational instability in 
which small low-momentum fluctuations in the polarisa- 
tion are amplified to macroscopic spin domains. In the 
miscible regime, a < 1, we invoke a counter-superflow 
instability J5l [16] by choosing oppositely directed flow 
vectors for the two field components. Above a critical 
relative velocity, momentum exchange between the two 
fluids causes the superflow to decay and spin domains to 
form. Subsequently, overpopulation at intermediate mo- 
menta is encountered for all values of the parameter a. 
Its microscopic origin is seen in the coloured distributions 
of the density imbalance in Fig. |2] which show strong 
fluctuations and the onset of domain formation. We esti- 
mate t\ 0:^ 10/cJi, with uji being the energy of the fastest 
growing unstable mode. Note, that the dynamics up to 
ti is different for different a and involves for example 
isotropisation in case of a < [M]. The further dynam- 
ical evolution of overpopulated states involves particle 
transport towards small momenta and energy transport 
towards high momenta [T7l[T8]. In Fig. [2] this is best ob- 
served by comparing the high momentum tails of the mo- 
mentum distributions at ti and ^ntpp • Particle transport 
towards low momenta eventually fills the zero-mode (not 
seen), a process most important in the miscible regime 
a = 0.8. The key insight gained from the spectra con- 
cerns the development of infrared (IR) power laws in the 
momentum distribution n{k) ~ at ^ntpp — lO^/cJi. 
These power laws depend on the external parameter a: 
( ~ 3.5 for a > 1, C 3 for a = 1 and ~ 4 for a < 1. 
Also, the pattern of imbalance fluctuations in the system 
depends on the interaction strength. 

We highlight the most striking observations at late 
times before we carry on developing a detailed under- 



standing of the relation between the dominant fluctua- 
tions and the observed power laws. In the immiscible 
regime the system consists of few large domains with 
additional inclusions of small point-like domains. As 
compared to the initial time domains have considerably 
grown. A special situation is encountered at the transi- 
tion point {a = 1), where domain-like structures persist 
to extremely long times. This is remarkable since in this 
regime domains are not energetically favourable as com- 
pared to overlapping particle densities. For a = 0.8, 
imbalance fluctuations have decayed up to few small ar- 
eas of strong imbalance. Let us finally look at the largest 
computed time — 10^/cJi. For a > 1, two domains 
of equal size remain, which reflects the immiscibility in 
the ground state. The small-scale domain- fluctuations 
have considerably reduced in number. At the transition 
point, we observe the persistence of domain-like struc- 
tures which we attribute to a diverging time scale r of 
domain decay as a ^ 1 from below [14]. For a < 1, 
small long-lived imbalance fluctuations remain, whereas 
the background tends to become very smooth, S*^ ^ 0. 

Finally, we investigate the microscopic origin of the 
scaling found in the momentum distributions. The bi- 
modal power laws in these distributions are signatures 
of the system having approached a non-thermal fixed 
point O |5] E]. At long times they become more pro- 
nounced in all three cases, signalling critical slowing 
down of the dynamic evolution. In Fig. |3| the result of 
a decomposition of the momentum distribution accord- 
ing to spin and fluid degrees of freedom at t = tNXFP is 
presented [14 . We show spin-spin correlation functions 
as well as correlations in the incompressible velocity as a 
function of momentum. The crucial point is that they ex- 
plain the scaling properties of the total momentum distri- 
bution and allow for an interpretation in terms of specific 
defect configurations, as shown in Fig.[2j This completes 
our understanding of transient non-equilibrium order in 
a two-component ultracold Bose gas. 

We remark that, for all a, compressible and quantum 
pressure excitations, which are not shown, dominate the 
spectrum in the ultraviolet (UV) regime and follow a 
thermal Rayleigh- Jeans distribution n ~ while they 
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give a negligible contribution towards lower momenta. 

In the immiscible regime, see Fig. [s] (left), the main 
contribution to the spectrum in the IR region is pro- 
vided by the incompressible component n^, corresponding 
to flow orientations transverse to the direction in which 
the flow velocity varies. Although we are dealing with a 
multi-component gas this feature is similar to superfluid 
turbulent flow in a single-component Bose gas. Thereby 
the incompressible spectrum shows ~ scaling over 
approximately one decade which is generated by coherent 
vortical flows around topological defects [Mj [171 [IH] • 
The cores of these vortex-like structures can be seen in 
the spin imbalance, see Fig. [2j since they are filled with 
particles of the other species and thus are of the skyrmion 
type. They are created during the merging process of 
domains and, persisting due to their topological nature, 
give the main contribution to the incompressible compo- 
nent spectrum. However, spin excitations overtake in 
an intermediate momentum region, showing a scaling of 
rig ~ k~^. Looking more closely, one observes that the 
scaling behaviour terminates in the IR at a scale tt/Ld 
given by the mean domain size Ld, while the cut-off in 
the UV at tt/^s is set by the width of the domain walls, 
i.e., the spin healing length ^5 = ^(1/|1 — Since 
the two contributions and rii are of comparable mag- 
nitude within an intermediate momentum range the sum 
of all contributions, giving the full spectrum n(/c), ap- 
pears to follow the power law n ~ k~^-^ in the IR, as 
discussed above. 

The situation on the other side of the transition can 
also be clarified. When a falls below 1 domains are en- 
ergetically suppressed. Thus, vortices dominate the non- 
thermal fixed point and induce the characteristic scaling 
n ^ Hi ^ see Fig. [s] (right). Here, it is also visible 
in the scaling of relative phase fluctuations nf^, which 
is related to constant particle densities on scales larger 
than ^ as compared to the skyrmionic case above. 

The picture changes dramatically on the transition 
point, a = 1. As shown in Fig. [s] (middle), vortical flow 
is much less important in this regime. The momentum 
distribution is dominated by spin fluctuations scaling as 
n ^ nl ^ ~ [M]. We attribute this feature 
to the increased dimensionality of the vacuum manifold 
for a = 1 which removes the topological protection of 
vortices. The stability of domain walls at the transi- 



tion point takes over but is non-topological in nature. 
Instead, conservation laws restrict the decay of this par- 
ticular type of defect [19]. This argument can be made 
even more transparent by studying the critical dynamics 
at the transition from a < 1 towards a = 1, where two 
non-equilibrium ordered states meet each other, vortices 
and domains [TT. 

We conclude that it is possible to establish the no- 
tion of a dynamical phase transition between different 
far-from-equilibrium ordered states in a two-component 
Bose gas. Beyond that, the underlying concept of 
(non-)topological defects determining bulk features of 
correlation functions in far-from-equilibrium situations 
is very general. It is easily imaginable that multi- 
component field theories with more than two components 
show similar behavior to the one outlined here. We point 
out the possibility of a universal duality between decay- 
ing defects and inverse particle cascades. This requires 
the generation of (quasi-) topological configurations far 
from thermal equilibrium and their slow decay, going to- 
gether with an increase of coherence and defect separa- 
tion [20]. Under these conditions, an inverse particle cas- 
cade is generated, and the associated power-laws can be 
found from the scaling properties of the respective single 
defect. A variety of (quasi-)topological excitations are 
known to exist in multi-component fields [TJ [19] , exam- 
ples are monopoles in gauge fields [ST and exotic mag- 
nets [22^, as well as skyrmions in Bose-Einstein conden- 
sates [ini |23] and liquid crystals [24]. New interesting 
features that are readily accessible in experiment are ex- 
pected for ultracold spinor gases [25 -29 . The transition 
between different types of transient non-equilibrium or- 
der can be controlled by changing the symmetry prop- 
erties of the Hamiltonian and thus topology and local 
conservation laws of the system. This offers interesting 
prospects for far-from-equilibrium dynamical phase tran- 
sitions in very different areas of physics. 
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I. SUPPLEMENTARY MATERIAL 

In the supplementary material we provide details of 
the methodology as well as additional results supporting 
our approach and conclusions. 

A. Classical field equations 

Since the dynamics we are interested in exclusively af- 
fects the low-momentum, strongly populated field modes, 
we employ the so-called classical field method which 
yields, within numerical accuracy, exact results for the 
time-evolving observables [1, 2 . For this, initial field con- 
figurations 01,2 (^5 ^o) are sampled from Gaussian prob- 
ability distributions and then propagated according to 
the classical equations of motion. At the end of the time 
evolution correlation functions are obtained from ensem- 
ble averages over the set of sampled paths. The classical 
equations of motion derived from the Hamiltonian den- 
sity (IT]) of the interacting two-component Bose gas are 
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FIG. 4: Momentum distribution n{k), defined in Eq. ([3|, at 
different grid times for the case a = 1. The counterflow mode 
is indicated by the dashed vertical line. The inset figure shows 
a snapshot of a typical spin configuration , corresponding to 
the earliest depicted spectrum at grid time t = 500. t = 5000 
corresponds to ti discussed in Fig. [2] in the main text. 



idt4>i = -lv^4>i + (l{\4>i? + a\<p2\^)4>i, (2a) 

idt4>2 = -\^^<p2 + 9{\<p2? + a|0i|')</'2 . (2b) 

The momentum distribution which is shown in various 
graphs is defined as 

n{k,t) = JdQk{rj{Kt)^j{Kt)). (3) 

where J dQk denotes the angular average in two- 
dimensional momentum space. 

We have applied the following rescalings to obtain di- 
mensionless quantities in the system of coupled Gross- 
Pitaevskii equations ([2|, using the lattice constant 
of the computational grid: 0, mg g and 

t/(mal) t. Here n = (A^i + N2)/L'^ is the mean total 
particle density on a A/'g x Ng grid of linear size L = NgQs^ 
with Ns = 1024. Conversion to physical time scales 
appropriate for a specific experiment, with, e.g.^ ^^Rb 
atoms, can be achieved by noting that the healing length 
^ = (2m^n)~^/^ is given by ^ = 4as for our system. 
Giving an estimate, e.g., for a physical healing length 
^ 1 /im, our lattice unit corresponds to 0.25 /im 

which yields a time unit of mRbttg/^ — 10~^ s. Hence, the 
characteristic time scale is 1/oJi — 10~^...10~^s which 
yields a realistic time scale for the appearance of the non- 
thermal fixed point at ^ntpp — 0.1 ... 1 s. 

B. Overpopulation from dynamical instability 

In our scheme the exponential growth of low lying mo- 
mentum modes is the first step towards the approach 
of the non-thermal fixed point. The production of an 



overpopulated state via instabilities has been discussed 
m Fig. [2] of the main text. Here, we review the initial 
dynamics in more detail for the case of a = 1, see Fig. [4] 
The counter-superflow instability leads to a growth of 
certain unstable modes, acting as a source for low mo- 
menta and depleting the counterfiow mode. This initiates 
isotropisation and a typical cascading process in momen- 
tum space, causing the transport of particles to the lower 
momentum modes while their energy is carried away by 
a few particles moving to the higher momentum modes. 
Fig. |4] shows this evolution until the earliest time shown 
in Fig. [2] in the main text when universal features start 
to dominate the dynamic evolution. 

C. Spin- fluid representation 

Although we always simulate the full set of equations 
([2| it is instructive to study the degrees of freedom 
which describe the relative evolution of the two com- 
ponents. Writing the fields in the polar representation 
4'i = ^/Pi^^Pi^^i) these are given by the local phase 
difference Or = Oi — O2 and by the local density dif- 
ference pi — p2- The Pauli matrices a'^ serve to define 
the Schwinger representation = (j)lafj(j)j (sum over re- 
peated indices implied) of angular momentum operators. 
This results in a three-component vector of (pseudo-)spin 
densities S^, a G {x^y^z} 

5" = 2v^cos(^,), (4a) 
Sy = -2^^ sin (Or), (4b) 
S'=pi-P2. (4c) 

where the modulus corresponds to the total density |S| = 
pi -\- p2 = pT' For convenience, we apply the redefinition 
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S'" PtS"" such that |S| = 1. Then, the Hamiltonian 
density H of Eq. ([T]) can be rewritten as 



Pt 



2p: 



+ 



9Pt 



gpT 



(i-«) [i-m'] , 



(5) 



which shows the influence of the relative degrees of free- 
dom and their coupling to the global ones [3 . The quan- 
tity Jt = PiV(^i + P2V(/?2 is the conserved total parti- 
cle current associated with a global U{1) symmetry of 
Eq. ([T]) associated with the global shift of the total phase 
Ot = Oi -\- $2. Thus Jt can not be expressed using just 
the spin densities but contains also the remaining fourth 
degree of freedom, the total phase Or- We obtain 



PtS' 



2 [1- (5^)2] 



{sy\/s'' -s'^vsy). (6) 



According to the representation ([5| of the Hamiltonian 
density the two-component Bose gas can be described by 
a spinor field S which is coupled to a hydrodynamic fiuid 
with density pr and (conserved) quasi particle current 
jx. For a fiuid at rest, i.e., pT = const and }t = 0, the 
spin system thereby assumes the form of a classical non- 
linear sigma model (NLaM) with a mass term ^p|./4(l — 
a) [{S^)'^ + (S*^)^], whereas in general a current jr 7^ 
leads to a highly non-trivial coupling between internal 
and hydrodynamic degrees of freedom. 



D. Hydrodynamic decomposition 

In order to gain further insight into the interplay be- 
tween the domain structure and other excitations we de- 
compose the energy spectrum according to Refs. [IHS]. 
The results of this decomposition at the non-thermal 
fixed point were discussed in Fig. [3] of the main text. 
In a density-phase representation, the kinetic part of the 
Hamiltonian density can be written in the following way: 



Pt 



V^^V^'^ + lwp. (7) 



The velocity field w is defined via the total particle cur- 
rent y/pT^ = jr, similar to the convenient choice for the 
one-component case. In this decomposition the first and 
the last terms are the quantum-pressure component and 
the classical hydrodynamic component, respectively. In 
contrast to a single-component fiuid the second term of 
Eq. ^ adds a pressure-like contribution to the kinetic 
energy which is produced by internal excitations only. In 
addition, the velocity field w and thus the correspond- 
ing part of the kinetic energy can be decomposed in a 
compressible and an incompressible part, w = Wc + w^, 
with V X Wc = and V • = 0, such that wave-like 
and vortical excitations appear in different parts of the 
decomposition, see Refs. [11 [5]. Based on Eq. ([7|), radial 
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FIG. 5: Top left: Integrated density imbalance or polarisation 
{\S^\) as a function of time for various a. Top right: Final 
values of (|*S'^|) as a function of a, showing a slow down of the 
imbalance decay when a ^ 1 from below. Bottom: Decay 
time t/tq as a function of a. 



energy spectra in momentum space corresponding to the 
discussed energy parts can be defined as 

Es{k) = ^ fdnk{\ws{k)\^), Se{q,c,i} (8a) 



(8b) 



Here we have introduced additional velocities Wq = 
V -^/pT and = ^/prVS^ /2 for the sake of a closed 
representation. Finally, the energy spectra can be con- 
verted to occupation number spectra defined in Eq. ([s]) 
by multiplication with a factor ns{k) = k~'^Es{k)^ 
S e {q, c, z, s} ^. 



E. Domain and vortex spectra 

In Fig. [3] of the main text we discussed the microscopic 
origin of characteristic power-laws in the momentum dis- 
tribution n{k) by decomposing the spectra in different 
components. This was successful since signatures of dif- 
ferent defects show up in specific components of the spec- 
tra. Here we make this correspondence explicit for the 
case of vortices and domain walls. 

Under the assumption of spatially constant to- 
tal particle density the spectrum Ug can be related 
to the Fourier transform of the correlation function 
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of the angle-averaged spin order parameter S{k) = 
JdQk (^»(-k)5'^(k)), 

k\s = ^ Jd^k (^[V5«*]^[V5°]) 

= p^e ldnk{S''{-^)S^{k)), (9) 

with denoting the Fourier transform, and therefore 
ns{k) = pTS{k)/2. Hence, in a momentum regime 
"k/Ld <C k <C 7r/^5 where the the scahng behaviour is 
dominated by a single domain wah, Ug ^ S fol- 
lows from an ansatz in terms of the Heaviside function 
5'^ = ±(1 - 2(9(x)), 5^ = 5^ = 0. This feature is simi- 
lar to the scaling induced by solitons in one-dimensional 
Bose gases, where the phase jump occurs in the bosonic 
field and induces a scaling uid ~ see Ref. [7 . 

Similarly, the appearance of quantised vortices induces 
scaling behaviour of the incompressible energy n^. On 
scales considerably smaller than the mean distance be- 
tween vortices, the effective velocity field |w^(r)| associ- 
ated with the rotational flow decays as 1/r with growing 
distance r from the nearest vortex core. As a result, 
the angle- averaged incompressible kinetic energy distri- 
bution ~ |w||/2 gives rise to the IR momentum spectrum 
ni{k) ~ whereas compressible excitations vanish. 



F. Domain decay near the transition point 



As discussed in the main text, the transition point 
a = 1 provides a particularly interesting effect. Despite 
their energetic instability domain walls constitute the rel- 
evant spatial pattern associated with the non-thermal 
fixed point with a scaling n(k) ~ k~^ . For further con- 
sideration, we show that this is due to a diverging do- 
main decay time r as a ^ 1 from below. To show this, 
we plot the time evolution of the integrated density im- 
balance = J dx^ \S^\/LP' as a function of time in 
Fig. [5] (top left) and at final time tq in Fig. [s] (top 
righty From this data, we can extract a decay time r for 
domain structures defined by r = to/|(|6'^|)/ — 
which is shown, on a double-logarithmic scale, in Fig. [5] 
(bottom). This definition of the decay time coincides 
with the first-order term in the Taylor expansion of an 
exponential decay but has the advantage to be applicable 
to non-exponential decays. One observes a divergence of 
r as a approaches the critical value a = 1 separating 
miscible and immiscible regimes. This opens interesting 
possibilities for studying critical exponents in a far-from- 
equilibrium phase transition [8]. 
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